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ABSTRACT 

We calculate the kurtosis of a large-scale density field which has undergone weakly 
non-linear gravitational evolution from arbitrary non-Gaussian initial conditions. 
It is well known that the weakly evolved skewness is equal to its initial value plus 
the term induced by gravity, which scales with the rms density fluctuation in pre¬ 
cisely the same way as for Gaussian initial conditions. As in the case of skewness, 
the evolved kurtosis is equal to its initial value plus the contribution induced by 
gravity. The scaling of this induced contribution, however, turns out to be qual¬ 
itatively different for Gaussian versus non-Gaussian initial conditions. Therefore, 
measurements of the kurtosis can serve as a powerful discriminating test between 
the hypotheses of Gaussian and non-Gaussian nature of primordial density fluctu¬ 
ations. 
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1. Introduction 


Redshift surveys of optically selected, and IRAS-selected galaxies have been re¬ 
cently used to estimate the distribution of galaxy counts on several scales (CfA and 
SSRS catalogs: Maurogordato & Lachieze-Rey 1987, Alimi et al. 1990, Gaztaiiaga 
1992; 1.2Jy IRAS survey: Bouchet, Davis & Strauss 1992, Bouchet et al. 1993). The 
distribution turns out to be non-Gaussian, the departures from it being more pro¬ 
nounced for smaller scales. This fact does not necessarily imply that primordial den¬ 
sity fluctuations were not Gaussian because nonlinear gravitational evolution does 
not preserve the shape of the fluctuations. Therefore, observed non-Gaussianity 
of the large-scale galaxy distribution may be either of primordial origin or due to 
subsequent nonlinear gravitational clustering, or the result of both. 

Simple versions of inflation naturally produce Gaussian fluctuations from the 
quantum fluctuations of the inflaton field (Guth & Pi 1982; Hawking 1982; Starobin- 
sky 1982; Bardeen, Steinhardt & Turner 1983). On the other hand, the density field 
generated by topological defects (Vilenkin 1985; Vachaspati 1986; Hill, Schramm 
& Fry 1989; Turok 1989; Bouchet & Bennett 1990; Bennett & Rhie 1990, Albrecht 
& Stebbins 1992) as well as in some versions of inflation involving multiple scalar 
fields (Allen, Grinstein & Wise 1987; Kofman & Pogosyan 1988; Salopek, Bond 
& Bardeen 1989) can be characterized as non-Gaussian. The comparison of pre¬ 
dicted distributions of counts with existing observations, in either of the Gaussian 
or non-Gaussian cases, may serve not only as the confirmation of the gravitational 
instability theory, but possibly as a strong discriminating test between the two 
classes of models as well. It is, therefore, of great interest for cosmology to compute 
the predicted statistics of counts in both cases of initial conditions. 

Numerous studies have been undertaken to calculate the gravity induced depar¬ 
tures from an initial Gaussian distribution. Analytical efforts have been mostly suc¬ 
cessful in the weakly non-linear regime, where perturbation theory can be applied. 
The higher-order reduced moments (combinations of moments measuring departures 
from Gaussianity) of the evolved distribution have been calculated (Peebles 1980; 
Fry 1984; Goroff et al. 1986; Bouchet et al. 1992; Juszkiewicz, Bouchet & Colombi 
1993; Bernardeau 1994a,b; Catelan & Moscardini 1994; tokas et al. 1995) and re¬ 
cently the whole distribution function (Bernardeau 1992; Juszkiewicz et al. 1995; 
Bernardeau & Kofman 1995). 

Given Gaussian initial conditions, the reduced moments, or cumulants, have 
been shown to obey the specific ‘clustering hierarchy’ (Fry 1984). In particular, 
skewness and kurtosis are expected to scale linearly and quadratically, respectively, 
with the square root of the variance of the density field, a. (We use the stan¬ 
dard statistical definition of skewness, s, and kurtosis, k, namely s = (5 3 )/a 3 , and 
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k = (S 4 )/cr 4 — 3.) This is precisely what has been measured in redshift surveys. Fur¬ 
thermore, in the APM angular galaxy survey, containing as much as ~ 1.3 million 
galaxies, Gaztaiiaga (1994) has recently found to exist a 2D (angular) analog of the 
3D hierarchy mentioned above. Bernardeau (1995), for Gaussian initial conditions, 
has indeed shown this angular hierarchy to be theoretically expected. 

A non-Gaussian field has non-zero intrinsic, initial cumulants, even in the limit 
a —> 0. Therefore, the observed linear scaling of skewness with a (Bouchet, Davis 
& Strauss 1992) has been used as a strong argument against non-Gaussian seeds 
of initial density perturbations (Silk & Juszkiewicz 1991; Coles & Frenk 1991). In 
particular, Silk & Juszkiewicz (1991) have argued that skewness is constant for 
the cosmic textures model, in conflict with the existing data. However, as pointed 
out by Coles et al. (1993, hereafter C93) and Luo & Schramm (1993, hereafter 
LS), a non-Gaussian initial perturbation must also undergo nonlinear gravitational 
evolution and the gravity-induced skewness may quickly dominate the initial term. 
Both groups have shown it to hold indeed true for some specific cases of non- 
Gaussian initial conditions: C93 using fully nonlinear N-body simulations, and LS 
in weakly nonlinear regime, applying perturbation theory. Fry & Scherrer (1994, 
hereafter FS) have calculated, in a simple and elegant way, the growth of skewness 
of the density field evolving weakly nonlinearly from arbitrary non-Gaussian initial 
conditions. FS have confirmed that in general the evolved skewness is equal to the 
initial, “increasingly unimportant” term plus the term which grows linearly with 
a. They have concluded that “observations of a linear dependence of the skewness 
on the rms density fluctuation therefore do not necessarily rule out initially non- 
Gaussian models”. 

This is true. However, is it also the case for the cumulants of order higher than 
skewness? In other words, does the nonlinear gravitational dynamics wash away 
all information about the type (Gaussian or non-Gaussian) of initial conditions? 
The first, and so far, single attempt to answer this question has been done by LS 
who have also calculated the weakly evolved kurtosis. Unfortunately, the paper is 
unclear on this subject. In the abstract, “the importance of measuring the kurtosis 
is stressed since it will be preserved through the weakly nonlinear gravitational 
epoch”. The formula (3.16) of LS for the evolved kurtosis, however, is a sensitive 
function of a\ This inspired us to come back to this problem in more detail. 

In the present paper, we compute the kurtosis of a density field, which has un¬ 
dergone weakly nonlinear evolution from arbitrary non-Gaussian initial conditions. 
More specifically, we calculate the perturbation theory lowest order relevant con¬ 
tribution to the evolved kurtosis in its most general form. The FS method of the 
real-space analysis of skewness can be applied equally well to the case of kurtosis. 
Therefore, we do so; we also follow mostly their notation. The organization of the 
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paper is as follows: in Section 2.1 we qualitatively discuss the effect of non-Gaussian 
initial conditions on the expected scaling of the cumulants of a density field with 
rms fluctuation. In Section 2.2 we derive the general formula for the evolved kur- 
tosis. In Section 2.3 we calculate the effect of focal biasing on this formula. In 
Section 3 we apply the formula to the scaled lognormal toy-model and give detailed 
predictions for the kurtosis - variance relationship in that case. Our conclusions 
are presented in Section 4. 

2. Weakly Nonlinear Kurtosis 

2.1. Qualitative Analysis 

We are interested in the statistics of a large-scale density field, described by 
the density contrast, 5 = [p(x, t) — (p)]/(p). For a Gaussian random process, all 
information is contained in the power spectrum, or its Fourier transform, the two- 
point correlation function, £ 2 - This, however, is no longer true for a non-Gaussian 
process. In terms of the p-point moments, the process contains non-vanishing p- 
point reduced functions, £ p , p > 2. By definition, (5) = 0. The next few moments 
are 


(<S(xi)<5(x 2 )) = 6(xi,x 2 ), 

(5(xi)5(x 2 )5(x 3 )) = £ 3 (xi,x 2 ,x 3 ), 

(<J(xi)5(x 2 )5(x 3 )<J(x 4 )) =£ 2 (xi,x 2 )6(x3,x 4 ) + 6(xi,x 3 )f 2 (x 2 ,x 4 ) + 

f 2 (x 2 ,x 3 )f 2 (xi,x 4 ) + f 4 (xi,x 2 ,x 3 ,x 4 ), 


(6(xi)6(x 2 )6(x 3 )5(x 4 )5(x 5 )) =^5 / (xi,x 2 ,x 3 ,x 4 ,x 5 ) = 

=6(xi,x 2 )f 3 (x 3 ,x 4 ,x 5 ) + cycl. (10 terms) + 
6(xi,x 2 ,x 3 ,x 4 ,x 5 ) . (1) 

In the above equation we have also defined an unreduced 5-point function, 

The reasons why we have done so will become clear in the next Section. For a 
homogeneous and isotropic random process, the p-point averaged functions depend 
only on relative distances. For example, when p = 2, £ 2 (x 4 ,x 2 ) = £ 2 (|x 4 — x 2 |), and 
so on. 

The cumulants of the one-point probability distribution function (PDF) are 
simply related to the p-point functions: the p-th cumulant, K'. p = £ p (0), where 
Cp(0) = £p(0,..., 0). Thoroughout this paper we will be interested in normalized 
cumulants, or JC p = JC p /a p , where 
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(2) 


a 2 = <5 2 ) 

is the variance of the one-point PDF. The third and the fourth normalized cumulants 
are called in the literature skewness, s, and kurtosis, k, respectively. For the fifth 
one, q , we coin here the name “pentosis”. They are related to the central moments 
by the following equations: 


_ (5=) 

d Q 5 

(3) 

<7 

_ (5 4 ) , 

A ^ 1 

(4) 

<T 

(5 5 ) 

= ^-/-lOs. 

(5) 


For weakly nonlinear density perturbations (a < 1) one can apply perturbation 
theory. Within the framework of perturbation theory one approximates the solution 
to gravitational dynamics equations for the density contrast as a series 

5 = 5(i) + 5(2) + 5(3) + ..., (6) 

where 5( p ) are found to be of order of [5(i)] p (Fry 1984; Goroff et al. 1986). The 
first order (linear) term grows exclusively by an overall scale factor, 5(i) ( x , t) = 
D(t) 5(x, ti), where 5(x, U) is the primordial fluctuation imprinted in the early Uni¬ 
verse at an initial time t % and D(ti) = 1. Hence, in the linear regime (a 1), 
the probability distribution preserves its initial shape. In particular, its normal¬ 
ized cumulants, e.g. the skewness and kurtosis, of equation (4) and (5), remain 
unchanged. When higher order terms in equation (6) become non-negligible, they 
cause the cumulants to evolve. Since the lowest order contribution is given by 

is”) = ([i (1 , + i ( 2, +.. .n = <[id )] p >+p <[i(i ) ] (, '- 1) i ( 2,> + ■... ( 7 > 

it leads to 


(5 2 ) = ([5 (1) ] 2 ) + 2 (5 (1) 5( 2) ) + 0([5 (1) ] 4 ), 
<^ 3 ) = <[^(i)] 3 ) + 3 <[5 (1) ] 2 5(2)) + ... , 

<5 4 ) = <[5 ( i)] 4 ) + 4 <[5 (1) ] 3 5(2)) W 0([5 (1) ] 6 ). 
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(9) 

( 10 ) 



If the initial fluctuation field is Gaussian the second terms on the right hand 
side of equations (10) and (8) vanish. These terms are odd products of the initial 
variables and their average for a multivariate Gaussian is zero. This ensures that 
kurtosis, equation (4), for a density field evolved from Gaussian initial conditions 
scales quadratically with a. However, for initially non-Gaussian fluctuations these 
terms cannot be expected to vanish in general. Therefore, in this case, we expect 
the lowest order contribution to the evolved kurtosis to scale linearly with the rms 
amplitude of the fluctuations. 

On the other hand, in the above sense, the skewness is not the best statistics to 
discriminate between Gaussianity and non-Gaussianity: its lowest order correction 
(eq. [9]) is of even order, so it is already non-vanishing for Gaussian initial con¬ 
ditions. Indeed, it is clear, from equations (8) and (9), as found by FS, that the 
evolved skewness must be equal to its initial value plus a term proportional to <r; 
only the coefficient of proportionality is sensitive to the type of initial conditions. 
Unfortunately, the value of this coefficient cannot be used as a test for their nature. 
While perturbation theory describes the evolution of the mass distribution, obser¬ 
vations probe the distribution of galaxies and there are many reasons for thinking 
that galaxies are biased tracers of mass. Fry & Gaztaiiaga (1994) have shown that 
if the galaxy density is a nonlinear but local function of the mass density then the 
galaxy density field possesses the same scaling hierarchy as the mass field (see also 
Juszkiewicz et al. 1995). Therefore, the agreement between the observed hierarchy 
of the cumulants of the galaxy field with the hierarchy predicted theoretically for 
the mass field is considered as a strong support for both assumptions of Gaussian 
initial conditions and local biasing. On the other hand, local biasing changes the 
hierarchical amplitudes. Any discrepancy between the derived from theory and 
the measured amplitudes can be, therefore, attributed entirely to biasing and it is 
impossible to deduce from them alone anything about initial conditions. 

In sum, as already noted by LS, kurtosis seems to be a better test of the type of 
initial conditions. This led us to investigate the evolution of kurtosis in some detail. 

2.2. Quantitative Formula: Derivation 

The first term on the right hand side of equation (10) is, using the definition (4), 

([5 {1) ] 4 ) = (3 + K i )af iy (11) 

Here, Ki denotes the initial (= first order) kurtosis. To calculate the next term, 
we need the second order contribution to the density contrast, 5( 2 ). Peebles (1980) 
derived 
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(12) 


^(2) — y e e ,a^,a + j ^,ap^-,ap 5 

where 

e(x) = «5 ( i)(x), and A(x) = J — (13) 

Hence, 

4 ([^(l)] 3 ^(2)) = y (O - 4 ( f3e ,a A ,a) + ^(e 3 A ia/3 A :Q/? ). (14) 

The first term in the equation above is simply 

20 . 5 . 20 , 5 , ir , 

y(e ) = — Qi <x (1) , (15) 

where q/ is by construction the initial, normalized, fifth central moment. It is 
related to the pentosis (normalized fifth cumulant) of the initial field, g?, by the 
equation 


Qi — Qi T 10 Sj, 

(see eq. [5]). The second and the third terms are 


(16) 


-4(e 3 e,„ A, a ) = -4 J (e 3 (x) e(x),„ c(x')) , ) — 


,a 


d 3 


X 


47r 


(e 4 (x)e(x 7 )), 


x — X' 


(17) 


and 


8 

7 


( e ^,a/?A,a/3) 


d 3 :r 7 d 3 :r 77 
47T 47T 


(e 3 (x)e(x 7 )e(x 77 )) 


1 1 




(18) 


Now we apply the identity Y. a /(x — x 7 ) = — V 7 , a /(x — x 7 ) to equations (17) 
and (18). The integral in equation (17) can be then integrated by parts, giving 

-4(e 3 e, a A ^ a ) =-q/a 5 {1) , (19) 

and equation (18) takes the form 
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( 20 ) 


8 

7 


( e 


d 3 x' d 3 x" 
An At: 


&'(0, 0,0,x / ,x / - 


1 1 



,ap 


In the above, we have used equation (1) and the property of translational invariance. 
To simplify calculations we will carry them on with the unreduced 5-point function 
, and will substitute the partition of £5 at the very end. 

Applying the identity 


V7 V7 1 3x a Xp 5 a p 4vr 

V a V^ n =- -5 - —5 a p d D (x ), 


x- 


( 21 ) 


where 5 d(x) denotes the Dirac 5 function, the integral in equation (20) can be 
rewritten as 


1 


ftVm+Cfe']. 


3 ~ (i) 

Here, the functional C is generally defined as 
f d 3 x' d 3 x" 3 (3 cos 2 0 — 1) 


( 22 ) 


C[/(x',x")] = J 


Air An 


x' 3 x " 3 


/(x',x"), 


cos 6 = 


x' • x" 


• (23) 


Let us introduce the normalized functional, C. defined by 


C[f] = C[f]/J( 0,0). 


(24) 


Equation (20) can then be finally expressed as 

= 5 (1 + (?[&']) a'ay. (25) 

Combining the equations (10), (11), (14), (15), (19), and (25), we obtain 

(d 4 ) = (3 + Ki) er 4 ^ + a' q / cr^, (26) 

where 

a ' = fi + (27) 

Similar calculations for a 2 = (5 2 ) yield (see eq. [18] of FS) 

a2 = a (i) + \ b Si tfi), (28) 
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where s t is the initial skewness and 


b = f + fcfe]- (29) 

From equations (4), (26) and (28), we obtain the final expression for the evolved 
kurtosis: 


n = Hi + L a + 0(a 2 ). (30) 

Here, 

L = a' q- - b (3 + Ki) Si , (31) 

with a 1 and b given by equations (27) and (29). 

Let us now recall that 6 7 = £5 + £ 2 ( 1 , 2) £ 3 ( 3 , 4, 5) + cycl. (see eq. [1]). Thus, we 
have 

^(O, 0, 0, x 7 , x 77 ) =6(0,0,0,x 7 ,x 77 ) + 

3 6(x')6(x // ) + 3 6(x // )6(x') + 

3 6(0) 6(0, X 7 , X 77 ) + 6(x' - X 77 ) 6(0), (32) 

with 6(x) = 6(|x|) = 6(0, • • •, x ). Let us calculate C[ 6 7 ] using the above partition 
of 6'- The first term is simply C[ 6]- The next two terms vanish by symmetry. The 
fourth one is 3C[6] a fi)- an d our calculation of the fifth one yields | 

By using equation (16) we can cast the formula for L, equation (31), in the 
following form: 


where 


L 


aqi + 


136 

~T 



Si 


(33) 


» = )( + b’lfel ■ 1111,1 * = fy + f ^Ks] • < 34 ) 

The formula (30) for the evolved kurtosis, with L given by equations (33)-(34), 
constitutes the main result of this Section. It shows that the weakly evolved kurtosis 
is not equal, as one might naively expect, to the initial kurtosis plus just the term 
proportional to the squared variance. Instead, as indicated by qualitative analysis 
in the previous Subsection, the term linear in variance is already present. 


2.3. Bias and Kurtosis 

Perturbation theory predicts the statistics of the mass fluctuations, while ob¬ 
servations probe the galaxy distribution. There are both theoretical arguments and 
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observational evidence for the idea that galaxies are biased tracers of the mass dis¬ 
tribution. Without discussing the general issue of biasing in detail, in the present 
Subsection we calculate the weakly evolved kurtosis of galaxy distribution, n g , in 
the local bias model (Fry & Gaztanaga 1994; see also Juszkiewicz et al. 1995). 

In this model, one assumes the galaxy density 5 g to be a nonlinear, but local 
function of the mass density, 

S g = f(5) . (35) 

The Taylor expansion for 5 g up to second order is 

S g = bS + \b 2 S 2 - \b 2 cr 2 , (36) 

where 

6 = /'(0), &2 = /"(0), <x 2 = <£ 2 >. (37) 

The last term on the right-hand side of equation (36) ensures that (S g ) = 0. The 
calculation of the second and the fourth central moment of S g , up to the lowest 
order nonlinear correction, yields 

= ( S g) = 6 2<j2 1 + ySjd + 0((7 2 ) , (38) 

and 

{S g ) = b 4 a 4 3 + Hi + ^L + 2 — {q[ — s*)j a + 0(cr 2 ) . (39) 

In deriving the above, we have made the substitutions (S 3 ) = Si<j 3 T 0(<r 4 ), ( 5 4 ) = 
(3TKj) cr 4 + L(j 5 + 0(cr 6 ) (eq. [30]), and ( 8 5 ) = g'cr 5 + 0(cr 6 ). The kurtosis in galaxy 
distribution can be therefore calculated, using definition (4) and equation (16), as 

Kg = + Lg a g + 0(a g ) , (40) 

where 

L g = — + 2 ^ Qi + (6 — Ki) Si , (41) 

and L is given by equations (33)-(34). Thus the local biasing preserves the scaling 
of the weakly evolved kurtosis. The coefficient L g is a function of the normalized 
cumulants of the initial mass fluctuations field and of the coefficients b and 62 . In 
principle, there are combinations of b and 62 for which L g vanishes. It is, however, 
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quite unlikely, because it requires a coincidental cancellation of the terms of different 
origin, namely the L/b term, arising from non-Gaussian initial conditions, by the 
62 term, arising from nonlinear biasing. 


3. Applications: Scaled Lognormal Model 

The simplest non-Gaussian models are purely local: a one-point PDF is non- 
Gaussian, but p-point correlations vanish, except at zero lag. Although physically 
not realistic, these models are very useful - due to their simplicity - to demonstrate 
the difference between the evolution of the kurtosis of initially Gaussian, and non- 
Gaussian density fields. This class of models has been investigated analytically by 
Fry & Scherrer (1994) and numerically by Messina et al. (1990) and Weinberg & 
Cole (1992). In the present paper we will restrict ourselves to the investigation of 
such models. 

In these models, C[£ 3 ] = C[^'} = 0. Therefore, to calculate the coefficient L in 
the general formula for the kurtosis (eq. 30), it is more convenient to use the form 
given by equation (31). The coefficients a' and b (see eq. [27] and [29]) are thus 
a' = 47/21 and b = 26/21. By using additionally (16), we obtain 

L= A [4? 9* + ( 392 - 26 «i) «<] • (42) 

As a model of one-point distribution function, we now use the (scaled) lognormal 
distribution. The PDF of the lognormal distribution, A(£, u), has the form: 


P(y) dy 



(log y - v) 2 
2£ 2 


dy 

y 


For the value of the parameter v we choose 


(43) 


v = -E 2 /2. 

The moments of such a distribution about the origin are 
/4 = TM"- 1 )/ 21 , t = exp [E 2 ]. 


(44) 


(45) 


The initial density contrast, e, is related to the random variable y by the equality 


e = c{y- 1). (46) 

The condition (44) ensures that (e) = 0. Indeed, (e) = c(y — 1) = c (/// — 1), 
which, from equation (45), is equal to zero. The coefficient c is uniquely determined 
by the condition 


11 



<t 2 } = <T, 2 


(47) 


which gives c = (r — 1) ^cq. Higher-order moments of e are given by 

(e n ) = c n (( 2 /-l) n ), (48) 

and can be easily expressed in terms of the moments of the variable y, equation (45). 
The resulting initial skewness, kurtosis, and fifth normalized cumulant are respec¬ 
tively: 


Si = (t- 1) 1/2 (t + 2), 

Ki = (r - 1) (r 3 + 3r 2 + 6r + 6) , 

q t = {t- 1) 3/2 (r 6 + 4r 5 + 10r 4 + 20r 3 + 30r 2 + 36r + 24) . (49) 

Here, r > 1 is a free parameter of the model (see the definition of r in eq. [45]). In 
this manner, we have constructed an initial field with a given initial rms fluctuation 
(T t <C 1. which can be arbitrarily far from Gaussian. This is different from the 
context in which the lognormal model has been introduced in cosmology: Coles & 
Jones (1991) have proposed this model as an approximation of the PDF evolution 
from Gaussian initial conditions. Indeed, for E —>■ 0, the distribution (43) tends 
to Gaussian. (Note that consequently the cumulants, equation [49], then tend to 
zero.) Therefore, to avoid possible misunderstanding, we refer to this model by the 
name Scaled Lognormal Model. 

The higher-order cumulants of our distribution are determined uniquely by the 
value of r, so anyone of them can serve as a parametrization of the field. In the 
following, we parametrize the initial field by the value of its kurtosis, . In Figure 1, 
we plot the coefficient L of equation (42), as a function of k*. (For simplicity, we 
investigate the non-bias case, so L g = L and n g = n.) If the density field evolved 
from Gaussian initial conditions, then L = 0 and the first nonvanishing term in 
the formula for the evolved kurtosis, equation (30), would be that proportional to 
a 2 . In the present paper we have not attempted to compute the proportionality 
constant of this term in the case of non-Gaussian initial field. Still, for reference, 
in Figure 1 we plot the value of this coefficient for Gaussian initial conditions, 
S 4 = 60712/1323 = 45.9. 

In Figure 2, we plot n/a 2 as a function of a. The dashed line represents the 
Gaussian initial conditions prediction, i.e. n/cr 2 = S 4 . The solid line is drawn 
according to the formula (30) for a non-Gaussian field, for Kj = 5. Here, S 4 has 
been accepted as the proportionality constant of the quadratic term. Finally, the 
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dotted line also corresponds to equation (30), but with the l La' term set to zero. 
From the plot one can see that it could be difficult to distinguish observationally 
between the dotted and the dashed curve, due to big errors in the small a part of the 
diagram (Bouchet et al. 1993; Gaztaiiaga 1994). On the contrary, the difference 
between the solid and the dashed curve is obvious. This illustrates the importance 
of the linear-in-er term in the formula for the evolved kurtosis from a non-Gaussian 
initial field. This term ensures that the kurtosis is very sensitive to the type of 
initial conditions. In particular, our toy model with /-q > 5 seems to be in clear 
conflict with observational data. 

It should be emphasized that our results hold firm in spite of the fact that 
we have not calculated the coefficient of the 0 (a 2 ) contribution to the evolved 
kurtosis. The departure of the value of this coefficient from the value we used, 
S 4 , can only shift the solid line in Figure 1 up or down, and in any case the solid 
and the dashed lines will remain distinctly different. (Even if we calculated the 
quadratic, 0 (a 2 ) term, the final controversy would remain about biasing, which 
can change the amplitude of this term for galaxy counts. Again, unless the biasing 
accidentally happens to cancel the linear, L term (see eq . [41]), the Gaussian, and 
the non-Gaussian initial conditions predictions will remain qualitatively different.) 

The lognormal distribution for the initial fluctuation field has been applied by 
Moscardini et al. (1993) to study, by means of N-body simulations, the evolution 
of the two-point correlation function from non-Gaussian initial conditions. In the 
present paper we are not trying to impose any constraints on the models investi¬ 
gated in that work, for a number of reasons. First, the lognormal distribution has 
been used by Moscardini et al. (1993) to model the gravitational potential fluctu¬ 
ation field, not the density one. Second, the initial field also had there lognormal 
correlations, absent in the toy model we investigated. And last but not least, we 
have performed our calculations for the case of an unsmoothed final field. However, 
in order to compare quantitatively the perturbation theory predictions with both 
N-body and observations, the effects of smoothing must be taken into account. 

The calculation of the effect of smoothing has proven to be a difficult task already 
for Gaussian fluctuations (Juszkiewicz, Bouchet & Colombi 1993; Padmanabhan & 
Subramanian 1993; Bernardeau 1994a,b; Catelan & Moscardini 1994; Lokas et al. 
1995). For a non-Gaussian initial field this is even more difficult, and we address 
this problem elsewhere (Catelan et al. 1995). Furthermore, the effect of the final 
smoothing of an initially non-Gaussian field is model dependent, since spatial aver¬ 
ages will then depend not only on the initial power spectrum, but on higher order 
reduced p-point correlation functions as well. Nevertheless, we think that smooth¬ 
ing cannot change the qualitative picture emerging from our calculations: being 
a linear transformation of the final field, this cannot add or remove terms in the 
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formula for the evolved kurtosis, equations (30)-(33). 


4. Conclusions 

In this paper, we have calculated the kurtosis of a density field which had un¬ 
dergone weakly nonlinear evolution from arbitrary non-Gaussian initial conditions. 
We have computed the perturbation theory lowest order relevant contribution to 
the evolved kurtosis. 

This contribution is in general linear in the rms amplitude of the density fluctua¬ 
tions, but it is absent for Gaussian initial conditions. Therefore, it induces different 
scaling of kurtosis with a from that for Gaussian fluctuations. We have shown this 
difference to be very distinct for the scaled lognormal toy-model. 

Galaxies are most likely biased tracers of mass. The local biasing preserves the 
structure of the Gaussian hierarchy of cumulants. The accordance with observed 
hierarchy in galaxy counts is therefore regarded as a strong support for both as¬ 
sumptions of Gaussian initial conditions and local biasing. On the other hand, the 
coefficients of the scaling (the amplitudes) of cumulants are substantially changed. 
Therefore, they cannot be used as a test of the nature of initial conditions. This 
fact additionally strengthens the importance of the different scaling of kurtosis. 

The subject deserves further study, the most important step being to take into 
account the effect of smoothing of the final field. However, as argued in Section 3, 
it is unlikely that smoothing can change qualitatively the picture outlined above. 
Still, even if it happened to be so, it would not be dynamics which wipes out the in¬ 
formation about the type of initial conditions! In the present paper we have shown 
that weakly nonlinear gravitational clustering preserves the information about ini¬ 
tial conditions, and that kurtosis is a good indicator of this fact. 
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Figure Captions 


Figure 1. The coefficient L, i.e. the coefficient of the iinear-in-er term in the formula 
for the evolved kurtosis (cf. equation 30) as a function of the initial kurtosis, for 
the Scaled Lognormal Model. For reference, the value of the coefficient of the term 
quadratic in a for Gaussian initial conditions (first nonvanishing one in this case), 
*S* 4 , is also marked. 

Figure 2. The evolved kurtosis over evolved variance, k/ct 2 , as a function of a 
for the Scaled Lognormal Model. The dashed line represents the Gaussian initial 
conditions prediction, i.e. n/a 2 = S 4 . The solid line is drawn according to the 
formula for a non-Gaussian field, equation (30) with L given by equation (42), for 
Ki = 5. Here, S 4 has been taken as the proportionality constant of the quadratic 
term. Finally, dotted line is also drawn from equation (30), but with the ‘Ler term 
set to zero. 
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